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I. Recent Progress Summary 


The research sponsored by NASA Langley Research Center Grant NAG- 1-7 
is addressing the reduced order model (ROM) problem in distributed para- 
meter systems (DPS) adaptive identification and control by focusing on 
the development and evaluation of an adaptive controller applicable to 
the active stabilization of a DPS, such as a large flexible spacecraft, 
given a necessarily reduced-order expansion approximation model structure. 

The adaptive control strategy chosen for investigation [1] was an outgrowth 

\ 

of earlier NASA-sponsored efforts to combine the modal expansion endemic 
to the large flexible spacecraft community [2] and self-tuning [3], which 
is the most widely accepted lumped-parameter system (LPS) adaptive control 
strategy. The annular momentum control device (AMCD) [4] [5] was used as 
the test example. 

The difficulty encountered in the initial work [6]- [8] supported by 
NAG-I-7 was the deleterious effect of control and observation spillover. 

Two interpretations of this problem, which is unavoidable in application of 
existing LPS adaptive control schemes to DPS expansion descriptions, have 
emerged. In [6] and [8] a discrete-time state model of a generic class 
of DPS 

V^(k) + V^(k) + f(k) (1) 

V R (k+l) = 0 M V N (k) + 0 R V R (k) + E r f(k) (2) 

y(k) = S* V k) + C R V k) (3) 

is formulated where the subscript N denotes the ROM quantities and R the 
residual model quantities. Assuming that the projection generating the 
ROM is done into the modal subspace yields 0 NR = 0^ = C). This state- 
space interpretation and attempts at its adaptive control resulted in the 
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following observations: (i) If there is no observation spillover, i.e. 

C R = (3, and no model error, i.e. 0^ R = (), then an indirect adaptive control 

strategy, e.g. [9], can be globally stable {8j. (ii) Unless C R = () and 

= 0. (1) is only a quasi autoregressive moving-average (AEMA) and not a 

true ARMA model [6], The alternate interpretation of the ROM effects 

presented in [7] is based on consideration of two types of identifier 

strategies. A DPS sufficiently accurately described by a very high 

dimension modal expansion 

N 

d(0,t) = I W. (t) 0 . (s) , (4) 

j”0 J 3 

where the 0j are the spatial mode 'shapes and their amplitude time 
histories, can be approximated by either a subset of the N modes of (4) 

M 

d (0,t) - l W. (c) « (s) ' (5) 

E j-0 J J 

A , A 

or by a selection of W, to minimize some measure of d-d with 

J 

M 

d(e,t) = l w (t) 0.(s). (6) 

j=o 3 3 

In neither (5) nor (6) must the M-KL modes used in the ROM be the first 
M+l in the full N+l mode expansion. The observations made in [7] were: 

(i) Since only d can be measured some signal processing is necessary to 
yield d^ in order to solve (5) for the nonphysical W. needed in the in- 
dividual modal dynamics identifiers and controllers. (ii) The fitting of 

A a 

d to d by selection of the VL may not yield adequately accurate values of 
the corresponding for stable modal identification and control. This 
latter fact was emphatically demonstrated in [7] by simulated use of (6) 
and resulting DPS instability. Observation spillover can be interpreted 

A 

as the source of difficulty in extracting d R from d or in having Wj 
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closely approximate if d is fit to d. Unfortunately the apparent use- 
fulness of proposed approaches to observation spillover reduction is 
severely limited [6J-[8]. 
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X I * Fut ur e P Ians 


Since observation spillover seems to defy complete removal , the question 
in applying available LPS adaptive control techniques to a ROM of a DPS 
appears to be their limits of tolerance to observation spillover* Un- 
fortunately, the LPS adaptive control literature contains very little 
generalizable insight into the reduced-order adaptive control problem* In 
fact the nonadaptive LPS reduced-order control problem is currently an open 
research question. 

Therefore in order to pursue the original objective of DPS adaptive 
controller development by judicious application of existing LPS adaptive 
control schemes, the thrust of the future research under NASA Grant NAG-I-7 
is being shifted to LPS reduced-order adaptive controller (ROAC) studies* 

The insight gained on this subset of the DPS adaptive control problem will 
hopefully prove transferable. Admittedly the ROM effects are slightly 
different. In LPS ROAC the problem is the lower order of the dynamic model 
selected; while in DPS ROAC the orders of the individual modal dynamics are 
assumed accurate but the number of inodes is insufficient. However the 


source of difficulty in both can be viewed as spillover* 

Consider a discrete- time LPS described by the transfer function 


H(z) 


Y(z) = b l + 

U(z) z-a 1 z ~ a 2 


(7) 


In a state description similar to (1)~(3) e would correspond to C R and due 
to the partial fraction expansion 0 D would equal 0 for the ROM 

N K 


?(k) = a y(k-l) + b u(k-l). 


( 8 ) 
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Assuming a stable control objective of 


G(z) = 1(f) = x < k > “ cx ( k ~D + dr(k-l) 


( 9 ) 


the control effort 

u(k-l) = f^rCk-l) + f 2 y(k-l) 

with 



( 10 ) 


(ID 


f 


2 


A 

c-a 

A 

b 


( 12 ) 


would adequately control (7) if e - 0 and b and a had converged to and 
respectively. 


Four generic possibilities will be tested for tolerance to nonzero e 
for this simplistic reduced-order model reference adaptive control (MRAC) 
task: (i) equation error based indirect adaptive control (El), (ii) equation 

error based direct adaptive control (ED) , (iii) output error based indirect 
adaptive control (OX) , and (iv) output error based direct adaptive control 
(0D) . The El approach uses equation error identification [10] of the 
plant parameters in a self-tuning [3] solution of the stated MRAC problem. 

The EO technique is based on an input matching [11] [12] solution to the 
MRAC problem. The OX approach uses HARF [13] as an output error plant 
parameter identifier in conjunction with the separation (or "certainty 
equivalent") assumption of self- tuning. The OD scheme is based on an 
output error identification interpretation [14] of the MRAC problem with 
a minor modification in the control law so the adaptive controller can 
be put in the form of the error system in [15] and proven convergent. 

Though these approaches are asymptotically equivalent for e - 0, 
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variations are anticipated in this ROAC test. The complete algorithms for 
each case are listed below: 

Plant: y(k) = (a^+a^yOc-l) - (a 1 a 2 )y(k-2) 

+ (b^+b 2 e) u (k-1) - (b^a 2 +eb 2 a^)u(k-2) (13) 


El: 


ED: 


01 : 


y(k) = a(k) y(k-l) + b(k)u(k-l) 

Kk+i) - a(k> + a lk-D froa-yoo ) 

l+yy .(h-i)+pu (k-i) 

b(bfl) - b(k) + . P . “(k-l>ly(k')-y(k)] 
1+yy (k-l)+pu (k-1) 

^-i (k) =— — 

b (k+1) 


(14) 

(15) 

(16) 
(17) 


L (lc) = — — ' L) 
b(k+l) 


(18) 


u(k) = ^(k) r(k) + f 2 (k) y(k) ‘ (19) 

p,r(k-l) [dr(k-l)+cy(k-l)-y(k) 1 

f (k) = f (k-1) + — r = _ — (20) 

(sgnb) (max|b | ) [1+p^r (k-l)+p 2 y 2 (k-l) ] 


a „ y 9 y(k-l) [dr(k-l)+cy(k-l)-y(k)J 

f 2 (k) = f 2 (k-l) + — 2 (21) 

(sgnb) (max|b|) [1+^r (k-l)+u 2 y z (k-l)] 

u(k) = f x (k)r(k) + f 2 (k) y(k) (22) 

y(k) = a(k) z(k-l) + b(k) u(k-l) (23) 


a (k+1) = a(k) ) ^ z ( k ~D MjOliM ± q(y(k-l) - z(k-l))] 
1 + pz 2 (k-l) + pu 2 (k-l) 


(24) 


£ (k+1) = b(k) + P u (k-D [y(k)-y(k) + q(y(k-l) - z(k-l))3 
1 + yz 2 (k-l) + pu 2 (k-l) 


z(k) - a(k+l) z(k-l) 4* b (k+1) u(k-l) 


(26) 
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(27) 


f x (k) - 


f 2 (k) = 


b (k+1) 

c-a(k+l) 


b (k+1) 

u(k) = f^k) r(k) + f 2 (k) y(lc) 
OD: x(lc) = c x(k-l) + d r(k-l) 


(28) 

(29) 

(30) 


f-^(k) = fj^k-1) + ii 1 r(k— 1> [x(k) - y(k) + q(x(k-l) - y(k-l))] (31) 
f 2 (k) = f 2 (k-l) + ji 2 y(k-l) [x(lc) - y(k) + q(x(k-l) - y(k-l))] (32) 
u(k) = f x (k) r(k) +.f 2 y(k) 

+ [1 + y 1 r 2 (lc-l)+n 2 y 2 (k-l)][x(k)-y(k)+q(x(k-l)-y(k-l))] 

' (33) 


Tha sensitivity of each scheme to £ will be compared via determination 
^ max ? ^ max „ 

of l (x(k)-y(k)) , k (x (k) -y (k) ) Z , £ u Z (lc), and k u (k) for various 
k k 

and b^. Clearly the relative damping of the modeled and unmodeled 
portions of the plant should prove important. The most basic ROM 

JL 

rule-of-thumb is the ability to neglect modes with settling times 10 
or less of the dominant modes. However, earlier research [6]-[8] 
indicates that the smallness of e will also be important. 
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ADAPTIVE CONTROL OF DISTRIBUTED IViKAXrTTER SYSTEMS: 
THE ULTIMATE REDUCED-OR t , v * pROBhrM 


by 


Mark J. 3aias 

Electrical & Systems Engineering Deportment 
Rensselaer Polytechnic Institute 
Troy, New York 12181 

A b s tract 

For many distributed parameter systems (DPS), such 
os highly mechanically flexible structures, it is essen- 
tial to provide stable on-board/on-line adaptive con- 
trol in the presence of poorly known system parameters. 
However, such a controller must be based on a reduced- 
order model of the 'DPS and spillover from the immodeled 
residuals con deteriorate the performance of the con- 
troller and, xn some cases, defeat the whole purpose of 
the adaptive control. 

In this paper, we investigate direct and indirect 
adaptive controllers based on reduced-order models of a 
DPS and point out the mechanisms whereby spillover con 
upset the stability of adaptive controllers. We present 
some conditions under which the adaptive -controllers re- 
main stable in closed-loop with the actual DPS and point 
out certain generic problems that must be overcome for 
successful operation of adaptive DPS control. 

1.0 Introduction 

A great abundance of adaptive control schemes ex- 
ists for linear, lumped parameter, small scale systems 
[l], whose system parameters are poorly known, and the 
application of adaptive control is being seriously con- 
sidered (and, in some cases, hotly debated) in a variety 
of areas ^ [2 J. The existing adaptive control methods can 
be divided into two (not unfriendly) camps: direct 
schemes, where the available control parameters are di- 
rectly adjusted (adapted) to improve the overall system 
performance, e.g. [3J-[63, or indirect (self-tuning) 
schemes, where the plant parameters are estimated by 
some reasonably fa3t system identification scheme and 
the control commands are generated from these parameter 
estimates as though they were the actual values, e.g. 

[7l [ 9 ]. The direct schemes are usually model reference 
adaptive , i.e. the controlled system ( pi aiitT is forced 
to behave like a model system which has the desired per- 
formance properties, such_ as transient response (pole 
locations). 

A 3 G, Kreisselmeier succinctly put it: "adaptive 
control trades plant uncertainty for uncertainty in the 
closed- loop system behavior”; the stability of a linear 
system which uses on adaptive controller is often in 
question because the closed-loop system will be highly 
nonlinear during adaption* Direct (or model reference) 
adaptive schemes have achieved a great amount of success 
in producing globally stable closed-loop behavior under 
certain technical restrictions on the plant [6]. Global 
stability of certain indirect schemes using adaptive es- 
timators has been shown in [8], [ 9 ] with fewer restric- 
tions on the plant replaced by the need for addition of . 
"sufficiently rich" external test signals to the control 
commands. Presently, all stable adaptive schemes are 
restricted to scalar (or at best, multi-input , single- 
output) systems; the stability for general multi variable 
systems remains an important unsolved problem in adap- 
tive control. Another fundamental restriction is that 
the plant must be of known finite dimension, and' this 
dimension must be small enough to meet computational 
constraints for the on-line use of the particular adap- 
tive scheme. 


C. Ricnard Johnson, Jr. 

Department of Electrical Engineering 
Virginia Polytechnic Institute & State University 
Blacksburg, Virginia 2^06 1 

The abundance of adaptive control senemea is over- 
whelming and an understanding of the interrelationships 
and structural commonality of these methods is desper- 
ately needed ; see [6], [lO]-[ll] for progress in that 
direction. 

This paper addre33es the problems that appear when"' 
even the most veil behaved (i.e. stable) adaptive 
schemes for lumped parameter systems (IPS) are applied 
to linear distributed parameter systems (DPS), i.e. sys- 
tems whose behavior must be modeled by partial differ- 
ential equations (PDE). It is natural to want to make 
maximum, use of the existing body of UPS adaptive con- 
trol theory in the new context of DPS and it is obvious 
that "something" will go wrong - according to a well 
known l oik theorem if it can go wrong, it will. We 
want t,*' 'point out some of the basic mechanisms through 
uhlt*h things can go wrong for adaptive control of DPS. 

Our viewpoint is shaped by experience with application 
of Dr:i control to large aerospace structures, e.g. £l2]- 
mid this paper is an expansion of some of the 
topics presented in [l?3 f however, what we have to say 
here is in the broader context of DPS where applications 
include such diverse topics as flutter suppression for 
aircraft [l6j and actively controlled civil engineering 
structures, e.g. tall buildings and long bridges, sub- 
ject to hign winds or earthquakes [lT]. 

Tim most fundamental, generic problem for adaptive 
D!^ control as that DPS are infinite dimensional and 
any inclement able adaptive scheme must be based on a 
reduced-order model (ROM) of the original system. There- 
fore, the adaptive controller order must always be 
smaller than that of the plant no matter how much on- 
line computational capability i3 available to implement 
iho coat 1 oiler. Tins violates the hypotheses of all the 
existing global stability results for direct or indirect 
adaptive control. In our estimation, this makes- adap- 
tive Olt* control the ultimate reduced-order problem and, 
until this problem is solved (or circumvented), the 
benefits of adaptive control will be denied to vhat may 
be its most needy custoner-D?5 applications. Prelimi- 
nary attempt^ at adaptive control for some specific DPS 
have Veen rade in [l8]-[2^]. 

2.0 DPS Deacrirtion 

The class of DPS considered here may be described 
by the l\.l 

v - Av + Bf; v(0) *= y (2.1) 

t u 

y « Cv 

where, for each positive tine t, the (possibly vector- 
valued) .system state v(t) is in H, an appropriate Hil- 
bert space vitn inner product (•,*) and corresponding 
norm |] * 1 1 . The operator A is a differential operator 
which ia time-invariant, has domain D(A) dense in H, 
and generate* a C Q se-aigroup U(t) on H. The control 

vector f is Itxl and the sensor output vector is Pxl and 
rank B=l! find rank C=P; this means that a finite number 

of control actuators(M) and sensors(?) axe used. 
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The semigroup U(t) is dissipative in the sense 

that 

| |u(t)| I < M 0 e" et , t > 0 (2.2) 

vhere ^ j> 1 and. e 0 (when £=0, we rasnure f!^-l). 

Thi3 means that the energy in v(t), i.e. j|v(t)|| 2 , is 
conserved or dissipated* It is known [2$j, [2(5] that 
if (2*2) does not hold with e > 0 trten it cannot be 
made so even with ideal 3tate feedback : 

f - Gv (2*3) 


case-Murphy 'e lav, again. 1 

The choice of subspace and type of projection is 

often a matter of clever design and insight into the 
specific DPS. However, there are some obvious candi- 
dates, e.g. Golerkin approximation using any appropriate 
spline basis [l?]. Another good candidate is the modal 
subacace* 

k s = sp {*,. V 

with 


because B has finite rank. Therefore, we will usually 
assume £ > 0 because most physical systems have some 
dissipative mechansims in them; however, £ muy be ex- 
tremely small (as it is in large aerospace structures). 
Equation (2.1) can represent a large class of DPS in- 
cluding mechanically flexible structures such as air- 
craft and spacecraft. 

3.0 Reduced-Order Modeling of DPS 

Even though. (2*3) niay be the Ideal feedback control 
lav, the full (infinite dimensional) state v is rarely 
available. Implement able controllers for most DPS must 
be based on on-board/on-line control computers which 
process the P sensor outputs y and produce the M 
control commands f. Hence, control design mu^t be 
based on a reduced-order model (ROM) of the DPS in 
( 2 . 1 ). 


a «k - Vv- (3 - 6) 

vhere ^ are the orthonomal mode shapes (eigenfunc- 
tions) of A corresponding to the eigenvalues X^. 

Such mode shapes are available when the operator A is 
normal and has compact resolvent; this is the case 
for most mechanically flexible structures modeled by 
(2.1). Note that, in many cases, the exact mode shapes 
may exist but not be known. When the modal subspace is 
used, the projection is assumed orthogonal and the model 
error terns and Ap, { are zero; however, the spillover 

terms depend on the actuator-sensor locations and the 
mode shapes and need not be zero. Other choices of sub- 
space and method of projection could be selected in an 
attempt to reduce the model error and spillover terns , 
e.g. [28]. 


R0M T s are produced by projecting the system (2.1) 
onto some ^ appropriate finite dimensional subs, pace H,j of 

H; such projections often are orthogonal but they need 1 
not be. Let v^=?v and v R =Qv where F is> the projec- 
tion of K onto H f j and Q is the projection of H onto the 
residual subspace (the subspace compJ ementary to in 
H ) and , from (2.1), we obtain: 

VVa + VR + 3 » f (3 - 1} 

- " 'V's + Vr + b r f (3 - ?) 

* = C N V H + C R V R (3 ‘ 3) 


4 . 0 Non adootive DPS Control 


When the parameters of the ROM ore completely known, 
a linear controller can be designed: 


f + JL , 


y + H 


12 


" u y * H i2 s 

H 2i y + K 22 S 


(U.l) 


where the controller gains are obtained from the 
ROM parameters (A^, C^) and the desired performance. 

Once the controller is designed, the effect of the re- 
siduals must be analyzed using (4.1) in closed-loop with 
(3.l)-(3.3). Energy bounds have been produced in [29] 
to predict the effect of spillover and model error on 
the stability and performance of the reduced-order non- 
adaptive control (4.i) with the DPS (2.1). 


where A^PAP, A^ R a PAQ, B R =*PB, etc. The terms 3 R f and 
C R v R are called control and observation sol Hover , re- 
spectively; the terms A,^ R v R and are called model 

error. The ROM is obtained from (3.1) and (3.3) with 

V ead c r =0: 

('ll* + V 0) 3Pv 0 (3 - U) 


Examples of controlled flexible structures have 
been produced [i3]» [30] where the spillover terms 
caused. instabilities in the residuals even though the 
ROM (3.4) in closed-loop with the controller (4.1) would 
have been stable. 

Digitally implemented controllers would be based 
on discrete- time versions of the DPS (2.1). On* such 
version can be obtained by using the uniform time step 
At with a constant control command f(k) over the inter- 
val (k^l)At£ t < kAt: 


The ROM state v„ and the residual state v tj fom the 

it K 

full state v» v^+v R » when the projections are ortho- 
gonal, 

I Ml 2 - l | v 3 || 2 - || v R || 2 ( 3 . 5 ) 

, All controllers designed on the ROM (3.4) must be 
evaluated in closed-loop with the full system ("5.1)- 
(3.3)* i.e. the effects of the residuals v, ; through 

model error and spillover must be considered. In some 
cases these effects are sufficiently I that they 
can be ignored and the DPS is controlled .juite well 
from a POM -based controller; most ti^es ctu'j'is not the 


{ 


v(k+l) « *v(k) * H f (k) 
y(k) = C v(k) 


vheie <? = U(At) and H * 


U(t) 3 dT. 


( 4 . 2 ) 


Other ver- 


sions could be obtained for example with non uniform time 
steps. Wheii the ROM procedure of Sec. 3.0 is used, the 
discrete-time versions of (3.l)-(3.3) become 


v.,(k+l) = V ; .{k) + * ;(R v R (k) + H.J f(k) (!».3) 

v R (k+:) = < * R . I v .;( k ) + $ R v R (k) + rtk) (i*. 1 *) 

y(k) => C„v,.(k) + C- v (k) (I*. 5 ) 

ti it r« n 
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5,0 Adaptive Control of DPS n’odel error tern are both zero, the quasi - APMA 

(5-2) is a true ARMA. for the DPS (U.2) and an/ stable 

The natural approach for adaptive control of DPS adaptive scheme based on this ARMA will be globally sta- 

aeema to be bi e vhen usei | in closed-loop with. the actual DPS (h.2). 


(a) choose a nice ROM; 

(b) use your favorite lumped parameter aaaptive 
control scheme; 

(o) design the adaptive controller as though the 
ROM were the actual system to be controlled; 

(d) analyze the closed-loop behavior of the adap- 
tive controller with the actual DPS. 

Often, step (d) is omitted. This can be disasterous as 
shown by the flexible struct lire example in [3l]. 

Even a globally* stable adaptive scheme may prove 
to be unstable when used on the full DPS instead of the 
ROM; this happens because of the effect of spillover and 
model error. The indirect adaptive control scheme of 
[83 uses an adaptive observer and would be globally 
stable on the ROM under certain technical restrictions. 
However, if it vere designed on the FOM (3.U) but used 
on the* full DPS (2.1), it would not necessarily be sta- 
ble; the mechanism whereby spillover enters the adaptive 
control scheme and spoils the stability proof is dis- 
played in [15]. This is not a failure of the method of 
[8]; it is a failure to satisfy the hypothesis of the 
stability result of PSl -everr DPS will fail tn satisfy 
this hypothesis unless the observation soil lover and the 
model error term are both zero 1 see Theorem in L15J). 

The model error tern is always zero when the exact mode 
shapes are used for the ROM but the observation spill- 
over need not be zero even then. Hence some means of 
compensation, such as an adaptive pre filter [l?3 or an 
adaptive version of the orthogonal filter [32], might 
be used to try to eliminate observation spillover, but 
this has not been fully studied. 

ABMA-Cettin 1 


The stable direct adaptive control schemes depend 
on an Auto-Regressive living Average (ARMA) representa- 
tion of the plant in discrete- time: 

y(k+N) » Z a r y(k+r-l)+ Z 6 f(k+r-l) (5*1) 

r*l r-1 r 

for some tf and appropriate matrices a^, 6^. What the 

ARMA says is that, after H time steps, the present out- 
put is only related to the post N outputs and inputs. 
Existence of an ARMA is directly related to the finite 
dimensionality of the plant (ll is usually that dimen- 
sion) and is obtained using the Cayl ey-Hami! ten theorem 
for matrices. For DPS, only* a '’quasi -ARIIA 1 ’ can exist. 
From the Appendix, we obtain the quasi -ARMA for the DPS 
(U.2) or ( U.3)-( !» .5) i 

f r H n f(k+r-l) (5.2) 


v(k+r-l) 

“r * r V 1IR ' ~ V'R 

and r r is defined in the Appendix. Since R(k)-0 when 

C =0 and A *=0, ve have the following result: 
n r 

Theorem: When the observation spillover ( and the 


H K 

y(k+tf) ** Z a y(k+r-l) + I 
r-1 r r=l 

+ R(k) 

where 


R(k) « C v (k+H) + Z A 
R R r=l r 


This result for direct schemes is analogous to bhe cue for 
indirect schemes in [15]. Both results require that 
observation spillover be eliminated somehow before it 
reaches the adaptive control logic. This is not easy 
to do in general i 


6.0 Generic Problems in Adaptive of DPS Control 

The most crucial problem of adaptive control of 
DPS is that the plant is infinite dimensional and, con- 
sequently, the adaptive controller must be based on a 
low-order node! *of the DPS in order to be implemented 
with an on-line/on-board computer. However, any con- 
troller based on a reduced-order model (ROM) must oper- 
ate in close- J oop with the actual DPS; it interacts not 
only with the ROM but also with the residual subsystem 
through tne spillover and model error terms. This con- 
tributes the following generic, difficulties for stable 
adaptive controllers: 

(l) interaction of the residual subsystem with the 
adaptive controller may negate the stabilizing 
properties of the controller unless observa- 
tion spillover can be completely removed be- 
fore it reaches the adaptive control logic; 

(?) the nonlinear nature of the adaptation mecha- 
nisms in these controllers substantially ag- 
gravates the residual interaction problem 
(note elimination of control spillover alone 
does not guarantee stability as it does in the 
nonad iptive case); 

(3) eo:"p«u nation for spillover (and model error) 
often 1 equires knowledge of the very plant pa- 
rameters that are poorly known-a vicious circle; 

(M the adaptation mechanism may shift the closed- 
loop frequencies around (i.e. the closed-loop 
system is time-varying as well as nonlinear); 
this may counteract any benefits from prefil- 
ter: ng to remove observation spillover. 

In adr.it ion to the above, the following also tend 
to spoil tne known stability results for adaptive con- 
trollers : 

(n) PPG 4n:*tr*il pioblems are often non-minimum 
phaii* due to noncOllocated actuators and sen- 
sors; this presents a stability problem for 
ui”cct scncr.es but not for indirect ones; 

(b) PI*" control must often be done with more than 
*'ne actuator and sensor; conversion of multi- 
variable systems to scalar systems via output 
feedback can be used to overcome the lack of 
a multivariable adaptive control but it intro- 
duces some new problems (i.e. it may destabi- 
lize the residuals) [l?]; 

(o ) indirect schemes need sufficient excitation 
from external test signals; however, these sig- 
nals may substantially excite the residual sub- 

Some robu.it or adaptive way to counteract observa- 
tion spillover must be found in order to make stable 
adaptive e- uirn! of DFG possible. Even then, observa- 
tion spillover will be present to some degree; conse- 
quently , global stability may be too muon to ask of 
these Development of spillover bounds, 

along the : ]-»■!, of those for nonadaptive controllers, 
nay lie able : gi /e nne idea of the regions of stabi- 

lity |V' m* •e.,i,f , ii operation of adaptive DF5 control; 
such 'le»pi *| ut ic 11** w in progress. 
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Appendix , Derivation of the Qnnsl-ARMA Ct.i n DrS 

In this appendix, ve will derive the quuni-APMA 
(5.2) for the DPS (MMt.s). 

Trcn (L.3), define h(k) 00 that 

v f| (k+l) = 4> ?J v {{ (k) + h(k) (A.i) 

Then, for any non-negative interne- L, 
t L r 

v J( (k+L) - 4 [( Vjj(k) + h(k*r-! : (A. 2) 

H»te that the sumaation in ( A.2 ) must be aero vhert L=0. 


R(K) = C p v R (k*S) + I Av R (k+r-l) 


A = r > - a c_ . 

r r HR r R 
This is the lets 1 red quasi-AKMA. 
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Abstract 

A recently developed strategy for adaptive sampled-data control 
of distributed parameter systems based on a plant modal expansion 
description and modal simultaneous identification and regulation 
algorithms is presented with frequent reference to the annular momentum 
control device (AMCD) test example. The requirement of observation 
spillover reduction, which is especially crucial to the proposed 
adaptive control strategy, is addressed. 
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NOTATION 


a, b modal differential equation parameters 

A number of actuators 

C number of modes controlled 

d distributed parameter system deflection 

d reduced order modal model deflection 

K 

approximate reduced order model deflection 
e identification error 

f actuator forces 

F modal forces 

M number of modes in approximate expansion 

N number of modes in accurate expansion 

s spatial variable 

S number of sensors 

t time 

T sample period 

W modal amplitudes 

a, 3 modal difference equation parameters 

Yj 6 modal controller parameters 

T AMCD spin rate 

0 ring particle angle with respect to reference particle 

X desired discrete characteristic equation coefficients 

P adaptive identifier step-size weights 



a 


actuator reference frame angles 
$ partial differential equation expansion spatial eigenvectors 

0} oscillatory modal amplitude time frequency 

ft set of distributed system particles 



Introduction 


The control of large flexible spacecraft has become an active 

1 

research and development topic, as demonstrated by a recent survey , 

Adaptive control of distributed parameter systems (DPS) is also an 

2-5 

emerging research concern * Since large ’ flexible spacecraft are 
acknowledged to be described by partial differential equations with 
uncertain, i.e, a priori indeterminable, parameters, such structures, 
requiring increasingly stringent shape and attitude regulation, are 
prime candidates for application of DPS adaptive control strategies , 

The objective of the present work is to provide a real-time 
simultaneous identification and control strategy applicable to DPS 
in general and large flexible spacecraft in particular. The real- 
time computation objective prompts the use of modal expansion descriptions 
of DPS to permit some parallel computation. The constantly activated 
adaptability of simultaneous identification and control is needed, 
for example, for adequate control of poorly behaved DPS, such as 
very lightly damped, large flexible spacecraft, during inadequately 

predictable plant parameter changes due to operating condition variations, 

4 

Toward this goal, this paper improves a previously proposed strategy 

t 

with further attention to a step-by-step adaptive controller development 

4 

procedure and its consequences* The original idea was to combine 
a truncated modal expansion description used in flexible spacecraft 
control"** with a simultaneous identification and control (also termed 
self-tuning^) adaption strategy to regulate the lumped parameter modal 


1 



amplitude descriptions. The procedure of the next section also 
mentions the possibility of direct rather than indirect modal control 

g 

parameter adaption, such as via model reference adaptive control . 

The outlined procedure also reacts further to the special problems 

9 

spillover creates for an adaptive implementation. This strategy 
is followed in .the third section to develop an adaptive, regulator of 
the linearised, out-of-plane deflection of a spinning annular momentum 
control device (AMCD) , a candidate for large flexible spacecraft 
control^’ Simulations of this application are presented in 
the fourth section. This example is also used in the fifth section 
to illustrate the difficult modal observation spillover reduction 
problem. An enlarged framework for adaptive control of distributed 
parameter systems, especially flexible spacecraft, is overviewed 
in the conclusion. ' Unfortunately, achievement of the grandiose 
objective stated at the start of this paragraph is only (possibly) 
n begun in this paper. 

A Distributed Parameter System Adaptive Control Strategy 


For DPS describable by 

CO 

d(s, t) = l W (t) $ (s), 
j=0 J J 


U) 


where d is a vector-valued function of spatial location s and time t, 
e.g. displacement for flexible spacecraft, <p ^ (s ) an orthogonal expansion 
basis of shape eigenfunctions, and W.(t) the amplitude of the j tb 

3/ A 

shape function at time t, an adaptive control strategy has been proposed 

Under the assumption that the amplitudes obey uncoupled, linear, ordinary 

differential equations of known order, but with unknown coefficients, e.g 


2 



( 2 ) 



paper. 

r 4 

The adaptive modal control strategy can be divided into two stages: 
one prior to system operation and the other on-line. The recommended 
steps of pre-activation analysis are: 

(i) Determine expansion basis <j>^ in (1). 

(ii) Select finite expansion upper limit to approximate (1). 

(iii) Specify sensor locations and relate distributed measurements of 
d(s, t) to the modal amplitudes W^(t) in (2) via reversal of (1). 

(iv) Determine actuator distribution and formulate effect on modal 


3 



forces F.(t) in (2) via (3), 

3 

(v) Establish modal control objectives* 

The recommended steps of real-time, adapting, sampled data control 
formation are: 

(vi) Apply previously calculated actuation forces and sense d(s, t) 
in (1) at sample instant. 

(vii) Process sensor data to estimate modal amplitudes fcL(t) in (2)* 

(viii) Select the modes requiring control. 

(ix) Process applied forces f(s, t) via (3) to determine achieved modal 
forces F^(t). 

(x) Improve the identification of the discretization of (2). 

(xi) Design modal controllers on-line with current parameter estimates 
to meet modal performance objectives. 

(xii) Convert desired modal control to actual actuator commands. 

(xiii) Repeat (vi)-(xii) at the next sample instant. 

The simultaneous identification and control strategy of steps (x) and 

(xi) could be replaced by a single step improving the controller parameters 

12-14 

by direct lumped parameter plant adaptive control' . The espoused 

indirect strategy is more intuitive and currently more flexible in 

terms of control objectives though more restrictive in terms of the 

13 16 

general need for plant parameter identifiability ’ as manifested 
in adequate model complexity and sufficiently exciting input requirements. 

As outlined the modal controllers rely only on past data to 
determine the present control action. This allows the considerable 
(despite the possible parallel execution of steps (x) and (xi)) 



computation to be done during the full sample period. The steps 

(vi)-(xii) differ slightly from their earlier description^ due to 

a fuller appreciation of the demands of steps (vii) and (ix), 

especially (vii) , which are discussed in later sections. Further, 

detailed comments on each stage of this step-by-step procedure 
4 

are available . The next section illustrates this strategy by 
application to adaptive regulation of the small, linearized, 
out-of-plane deflection of a large, flexible AMCD. 

AMCD Application 

Large momentum vectors resulting in the rotation of large space 
structures can be created smoothly by a solar-powered (and therefore 
effectively non-depletable) , dual momentum vector configuration of 
two counterrotating AMCDs magnetically attached to the space structure^ 
The AMCD components of these attitude control devices will probably be 
as large in diameter as possible for such structures and as small in 
cross section as possible in order to maximize their momentum/mass 
ratio* Therefore such AMCDs would behave like lariats with translation, 
rotation, and deformation modes of disturbance from their nominal planar 
spinning configuration. 

The step-by-step procedure of the preceding section will be 
followed in designing an adaptive modal controller of such an AMCD. 

(i) The boundary conditions of ring closure permit the use of a 
Fourier series to describe the linearized, out-of-spin-plane deformation 
of the AMCD 


5 



d(6, t) = l [W (t) cos(je) + W (t) sin( j8) ] , (4) 

j-o CJ SJ 

which is of the form of (1) with the angle 6, measured around the AMCD 
ring from a reference (0=0) particle to the location in question, 
as the single spatial variable s. This sinusoidal basis is also the 

17 

eigenvector basis for this out-of-plane motion of a homogeneous -ring . 
(ii) Assuming that the higher spatial frequency deformations will 
exhibit lower amplitudes, (4) can be truncated with arbitrary accuracy 
as 

N 

d(0, t) = £ [W (t) cos(j6) + W (t) sin(j0) ] . (5) 

j-0 SJ 

This limit N may permit accurate approximation of (4) but be an infeasible 
limit in- terms of controller computations. If a further reduction is 
necessary, ‘(5) and therefore (4) (and in effect (1)) will be roughly 
approximated by either 

M of [C>, N] 

d (0,.t) = l [w .(t) cos(jG) + W (t) sin(j0)j, (6) 

K Cj Sj 

M of [0,-Nj 

where £ signifies a summation over the index j where j is any 

3 

Mfl entries of the set {0, 1, 2, .... N}, or 
M of [0,N] . 

d(e, t) = l [W (t) cos ( j 0) + W (t) sin(j6) ] , (7) 

j C J S J 

- A 

where the are not necessarily the corresponding in (6) and (5) 
due to the selection of W to best fit d to d given 2M + 1 point 
measurements ,of d, e*g, as in the Galerkin approach^. 


6 



(iii) Assume S sensor measurements of ring particle deflections 
d(0 i9 t) at i 3 1, 2, , S can be processed simultaneously. These 

measurements can be decomposed into modal amplitudes W. by multiple 


concatenation of (6) 
for the lower modes, 


M of [0,N] 

(or (7)) with £ replaced 

3 

as 


M 

by l , 

j=o 


i. e. 


V e i*t) 


d R (6 2’ t) 

= 

: # 

♦ 

vv* 



1 cos(0^) ... 

1 cosC©^) ... 

1 cos(0g) ... 


cos(M0 ) sin(0^) 

cos(M0 2 ) sin(0 2 ) 

• ♦ 

* ♦ 

• • 

cos (M6g) si‘n(0g) 



• (8) 


With appropriate reindexing of the right side of (8) any M modal amplitudes 

composing d in (6) could be written in this matrix form. The 0. may 
R i 

vary from sample to sample, especially if the sensor (s) is not spinning 

with the AMCD. Note that (8) requires measurement of d not d. Obtaining 

K w 

d^ from d requires observation spillover removal, as will be addressed 
later . 

(iv) Assume A actuators are located in a reference frame fixed with 
respect to AMCD spin, e.g. the suggested'*'* magnetic "bearing" actuators 
attached to the spacecraft. From (3) where ft = { 0 [ 0 s [0, 2 tt) } an 
assumption of* point actuation located at converts the integrals to 
summations over the set of A actuators. Note that, due to an AMCD spin 
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rate of T radians /second relative to the actuator locations, the a. 

1 

must be converted to the ring particle reference frame via 


e i (t) = CT i ( t ) ~ rt > 


( 9 ) 


which assumes that the references 0=0 and o - 0 were aligned at 
t = 0. Similar to (8), for the C modes to be controlled 


1 

o 

/—N. 

rt 


1 

1 

i 

F c i<t) 


cos(Q^) 

cos( 0 2 ) ... 

f 

cos(0 A ) 

F cC^ C> 


cos(C9^) 

» 

cosCCS^) — 

cos(C 0 ) 

A 

F sl (t> 


sin( 0 1 ) 

sin( 0 2 ) ... 

• 

sin(0 A ) 

F sC (t > 


sin(C 0 ^) 

9 

9 

sin(C02) • « * 

m 

sin(C0 ) 

ri 


f( 0 l,t) 




( 10 ) 


can be formed. 

(v) Ring stabilization requires mode damping. Satisfactory modal 
damping can be provided by modal pole placement. 

(vi) If d(0, t) is provided by an effectively instantaneously scanning 
sensor then any S ring particles can be observed. If the sensors are 
incorporated with the actuators then (9) must be used since different 
ring particles will be sensed at each sample instant. One other 
possibility is more frequent sensor interrogation than actuator 
reactivation, allowing additional signal processing 'possibilities 
before control selection. 
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(vii) Solution for the [W ] vector in (8) can be achieved by 
pseudoinversion^ or by a DFT^ if the 6^ are equally spaced. Note 
that if the 0^ are equally spaced such that 0_^ = 2iri/S and S = 2M, 
then. M0_^ = iri and the rightmost column of the S * (2M + 1) matrix 
in (8) equals zero. Therefore the column of sin(M0_^) entries should 
be removed to retain invertibility. In such a case the Mth mode 
sine component is unobservable. 

(viii) If C < M, the current strategies are to control those C 
modes either with the lowest spatial frequencies or with the greatest 
modal amplitudes. The assumption leading to (5) will tend to equate 
these two classes. 

(ix) If A < 2C + -1 in (10) then due to the least squares solution 
implemented in step (xii) of the last sample instant the desired 
most likely will not have been achieved and (10) must be calculated 
to determine the forces actually reaching each mode needed in the 


next step. 

(x) For a large, lightweight AMCD (2) will be second order and 
essentially undamped 
2 

[W.(t)] + u, 2 (F, t) W (t) = (tj-) F (t), (11) 

dt 3 3 3 j 3 

where the M denote the effective modal masses and the T and t arguments 

of (o are intended to evoke the slowly time-varying character of the modal 

amplitude time frequency due to such operating conditions as spin rate 

and temperature. Note that the magnitude of the oscillatory initial 


condition response of (11) is inversely proportional to ai.. 


The notation 
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of ( 11 ) (and subsequently ( 12 )-( 20 )) is intended to encompass both 

W . and W .. Note that both W . and W . have the same to , i.e. 
c 3 sj cj S 3 3 

oj . = w .. Assuming uniform sample intervals of T seconds and 
C J S3 

constant modal forces over the sample period, 

WjOO = 0jl W j (k-1) + a. 2 V. (k-2) + g j;L F j (k-1) + $ j2 F (k-2) (12) 

is an exact discretized predictor of the modal amplitude where 

cx^ = 2 cos (u^T), = -1 and = P ^ 2 = (1 - cos W..T) /(M^ 2 ) . 

The-addition of damping in (11) will effect the definitions of o 

and g but not the form of (12) . Note that constant actuator forces 

f over the sample interval will not generate constant modal forces 

due to AMCD rotation. This can be combatted by actuator force 

windowing. If the modal forces vary over the sample period, (12) 

becomes an approximation only as accurate as the degree of constancy 

of Fj (t) over (k-2)T < t < (k-l)T and (k-l)T < t < kT. If the 

structure of ( 12 ) is used for an adaptive identifier or to structure 

a direct adaptive controller then uncertainty in ok, M , and the 

neglected damping coefficient can be accomodated. The anticipated 

problem in AMCD control is uncertainty in . The second-order (12) 

can be identified by either of two broad classes of recursive parameter 
?0 

estimators termed - prediction error and pseudo linear regression and 

21 22 

represented by equation error and output error identifiers, 
respectively. An equation error formulated identifier for (12) of the 
form 
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“jl (k> 


V k) 


V k) 


V k > 



a (k-1) 
a^ 2 (k-l) 

6 2 (k-l) 


e (k-1) 


1+ l []i (k-l)W?(k-l-l)+p,.Oc-l)F?(k-i-l)] 


where 


U (fc-l)W (fc-2) 

U (k-l)tf (k-3) 
3 o 3 
p (k-1) F (k-2) 

p, 2 (k-l)F.(k-3) 


(13) 


e. (k-1) *» W. (k-1) 
3 3 


l [a . . (k-l)W . (k-i-1) + g . . (k-1) F . (k-i-1) ] (14) 

i=l 31 2 31 3 


and 


0 < y . . GO < y (k-1) 

J 1 J 1 


<2 and 0 < p . . (k) < p . . (k-1) < 2 
Ji Ji 


V i, k, (15) 


requires exact measurements of sufficiently rich F_. and W^. for consistent 
identification. 

(xi) Feedback regulation structures of second order dynamic output 
feedback or equivalently (in the absence of unmeasurable inputs) 
reconfigured state variable feedback achieved via 


(k) = Yj i (k)W^(k-l) + y_. 2 (k)W^ (k-2) + 6 (k) (k-1) + 6j 2 (k)F^(k-2) (16) 

will cause the modal plant-^controller characteristic equation to converge 

to jz** + X + A. 0 z + X if che controller parameters are chosen 

3l J2 J]3 j4 
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via 


8 jl (k) - ~ « j]L Ck) (17) 

Yj X (k) = [6^ 1 (k)a^ 2 (k) - A^ + (S^OOct^U) - a_. 2 (k) - A^) 

* 

*(a (k) - a (k)0 (k)/B (k)) + A 8 (k)/g (k)]/ 

Jl 3 2 Jl J 2 3 4 Jl j2 

[^(kJa^Ck) + 3^2^ " ot^ 2 (k)$^ 1 (k)/B_. 2 (k)] (18) 

5 j2 (k) = ywyk) - A j2 - 3 jl (k)^ J1 Ck) - a 2 (k) (19) 

Y j2 (k) = (6 j2 (k)a j2 (k) - A j4 )/6 j2 (k) (20) 

^ A A 

and a + a and g g. Note that the strictly causal formation of a and 
8 in (13)-(14) used in (17)-(20) to parameterize (16) permits the 
assumed predictive control formation. Furthermore, the form of (16) 
purposely avoids velocity measurements, which are expected to be 
difficult to sense for large AMCDs. 

(xii) Least-squares solution of (10) for the f(*,t) Is required if 
A < 2C + 1. 


Simulations 

Successful simulations of this AMCD example have been reported 
4 11 23 24 

elsewhere J * 5 Consider here' an AMCD described by (5) with 

N = 4 and the w_. in (11) equal to 0, 2.62, 5.24, 7.85, and 10.47 for 
j=0, 1, 2, 3, 4 respectively. For an AMCD spin rate of T = 60 
degrees/second, a sample period of T = 0. 1 seconds, and the same 
damped pole placement objective for each mode of A^ = -1.684, 
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A^ = 1.165, A^ = -0.402, and A^ = 0.0558, the strategy of the preceding 
section was applied in the following cases: N = M = C = 4, 

N = M= 4>C=2, and N = 4>M=C=2. For each situation the 
AMCD was given an initial deformation composed of unity modal amplitudes. 
The adaptive controller was applied, with initial modal frequency 
estimates . = {4, 8, 12, 16} converted to initial parameter 

estimates in (12) via the discretization formulas, in an attempt to 
stabilize the ring deflections to zero. 

The AMCD simulation consisted of (10) with C replaced by N, which 
converted the applied forces f to modal forces F which were assumed 
constant over the sample interval, the use of the appropriate F in 
(12) to update the N modal amplitudes W, and the formation of the 
deflection d via (5) for all 6. The controller began with measurements 
of d for use in (8) , which were supplied by solution of (5) at the 
appropriate sensor locations 9. The sensor locations were assumed 
to be rotating with the AMCD. A special reflective mark, on the AMCD 
and a centrally located, scanning laser detection system can be 
hypothesized as providing such ring particle deflections. The 

measured d were used for the d in (8) with M * 4 or 2 as required. 

K 

Solution of (8) provided the W used in (13) and (21). The improved 
parameter estimates provided by (13) with all p = p = 1 were used in 
(22)-(25) to parameterize (21). When M > C, only the F for the 
modes of lower spatial frequency were calculated. The past F in (21) 
were provided by the previous results of (21). Step (ix) could be bypassed 
since A was always chosen equal to or greater than2C + l. Inversion solution 
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of (10) provided the applied forces f to the AMCD simulation* 

Figure 1 with N * M - C = 4 illustrates the anticipated 

effectiveness of the adaptive strategy in recovering from inaccurate 

estimates of the co_. and successfully stabilizing the initial AMCD 

deflection, eventually to zero displacement. The asterisks on the 

displacement curves mark the measurement points which are fixed 

in the ring particle reference frame. The arrowheads along the 

spatial coordinate 0 axis show inertially fixed actuator locations 

with the length of the arrowshaft proportional to the applied force 

according to the right-hand scale. The plots are drawn in the ring 

particle reference frame so the actuator locations appear to regress 

for the progressing ring. The success in figure 1 is not universal 

for N - M - C. With only the desired pole locations changed to 

values nearer the unit circle by A^ - -2.314, A^ = 2.429, A^ - -1*176, 

A^ - 0.2, the proposed adaptive regulator fails via an often neglected 

stall mechanism. For example, if the parameters in (13) when used 

2 

in (17) -(20) lead to an unstable modal controller, W . will become 

J 

so large in the denominator in (13) that the parameter estimates 

are only insignificantly corrected. If this condition persists 

long enough, all claim of linearity can be abandoned in application, 

thereby negating claims of eventual identifier convergence. This 

stall mechanism, a characteristic of simultaneous identification 

and control not peculiar to just the proposed adaptive DPS control 

strategy, occurs in the cited example with the W ^ (5 seconds) = 

4 

4.8 x IQ and growing. This failure can be attributed to the 
smaller stability margin of higher modal frequency estimates for 
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low frequency, light: damping objectives. Different initial parameter 
estimates nearer the actual modal frequencies can avoid the stall of 
this particular example. Also the sufficient excitation requirements 
of simultaneous identification and control 7 are almost surely not 
met in this regulation example, thereby encouraging an expectation 
of failure. However, this pathological case does emphasize a need 
for an understanding (currently nonexistent even in lumped parameter 
system adaptive control) of closed-loop singularity migration during 
simultaneous identification and control. 

Figure 2 illustrates the boundedness of the AMCD deflection 

achieved in adapting from inaccurate ok prespecification with 

<v N = M=4>C-2. The visible higher spatial frequency remnants 

at t = 15 seconds are due to the, choice to control only the lower 

spatial frequency components of d, which do decay to zero. Asymptotic 

regulation of the full deflection is not achieved despite convergent 

parameter identification in the satisfactorily complex model due 

to control spillover arising from use of a restricted complexity 
9 

controller . For lumped-parameter systems this limited success 

25 

of adaptive control when N = H > C has also been documented elsewhere 
For A > 2C + 1, i.e. an oversufficient number of actuators for 
solution of (10), pseudoinverse solution of (10) provides a minimum 
control f energy solution, which can be expected to reduce the 
deleterious control spillover. Figure 3 with A = 9 for C = 2 
versus figure 2 where A = 5 for C - 2 clearly demonstrates this 
effect. 
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The possible failure of the proposed adaptive control strategy 

when N > M is documented in figure 4 where N=4>M=C = 2. The 

identifier of (13) is unsuccessful in converging on the actual a and 

6 values due to the use of d and not d in (8). Since d is used 

R 

in (8) and a least squares solution is used to determine the W, 
then actually (7) and not (6) is applicable. Figure 4, therefore, 
shows the failure of attempting to use a modal control strategy when 
approximately matching the full AMCD behavior with a reduced number 
of modes ♦ For this example, even if the a and 8 were successfully 

/V 

identified, using the W. and not the W. in (16) leads to instability. 

3 3 

Clearly a mixture of full behavior estimation and reduced-order 

control strategies, though commonly pursued in practice, is only 

valid if the modes omitted from the model contribute negligibly to 

the total behavior. Due to the control spillover reaching the 

unmodelled modes, even in the event of no initial energy in these 

radoes this negligibility can not be assumed. The difficulty adaptive 

control experiences due to N > M has also been noted in the lumped 

26 

parameter system case * The next section considers a signal 
processing strategy to combat this possibly- fatal problem by 
filtering the d of (5) to provide the d^ of (6) to (8), 

Observation Spillover Reduction 

Inexact sampled W^. are provided to the modal identifiers due 
to two sources of error: (a) aliasing, both in time and spatial 

frequency domains, due to discrete measurments of d(s, t) and 
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(b) reduced-order modeling inaccuracy due to N > M. However, 

based on the characteristic of flexible spacecraft that higher modal 

frequencies, e.g. the higher frequency for the AMCD in (4) as 

j becomes large, have correspondingly higher modal amplitude time 

frequencies, i.e. to. > to. for i > j, a strategy for extracting d 

1 j K 

24 

in (6) from measurements of d in (5) has been postulated . This 

approach assumes (a) satisfactory time-based sampling to avoid 

aliasing of modal frequencies up to w , i. e. T < — , (b) a sensor 

system capable of interrogating any ring particle at any sample 

# 

instant, (c) band-limiting knowledge of the M sought ok, e.g. for 

the lowest M consecutive ock specification of a frequency comfortably 

between a) and to , and (d) frequency-limited spectra for the 

F,. leading to separable W spectra at desired cutoff points. 

For a lowest M frequency approximation, implicit in (8), the 

following strategy appears reasonable- Assuming equally-spaced 

sensor measurements, in both space and time, and S - 2M, as in 

point (vii) of the preceding section, a DFT could be used to solve 

(8); however d and not d must be available. If the same ring particle 
K 

can be measured for deflection .at successive sample instants despite 

AMCD rotation then the sequence {d(0^, kT) } over k for a particular 

i, which is proportional from (5) to a fixed weight sum of the W^(kT), 

can be low-pass filtered (LPF) between and a> ^ removing high 

time, and therefore also spatial, frequency components leaving 

(d (8., kT)}. For each time t composition of (8) is now possible. 

K X 

The necessary assumptions preceding this particular strategy 
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description are quite restrictive. Due to the true form of (1) and 
the possibly tremendous magnitude of N and therefore (a) is 
mathematically impossible and only marginally practical. The limitation 
of (b) is that the sensors can not be colocated with the non-rotating 
actuators. Even if AMCD spin could be accomodated for DFT purposes 
the same ring particle could not be sensed at sequential sample 
instants disallowing the benefit of LPF. Conversely, always measuring 
the same ring particles could, in certain cases, lengthen the lag 
time before reaction of the control system to slowly propagating 
localized disturbances. Assumption (c) may be reasonable for 
discrete sinusoidal spectra but since (11) is forced by a nonlinearly 
generated, nonstationary (during adaption) signal in violation of 
(d) , the region of non-overlap between and becomes so small 

(if not nonexistent) as to require highly refined a priori knowledge. 

LPF phase distortion must also be assumed negligible. Clearly some 
higher spatial frequencies will be rejected by the LPF, but due to 
dissatisfaction of the assumptions pointing to unsatisfactory spectrum 
separation, d„ will be inexactly obtained. These reservations are 
even more severe for more complex filtering schemes, such as the 
effective comb filtering suggested earlier^ to be achieved via 
phase locked loops. 

One alternative^ satisfying the (c) and (d) requirements for 
filtering d_ from d is identification of the free AMCD response. 

However, this approach does not meet the stated simultaneous 
identification and control objective . Using a modified gain scheduling 
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concept to provide a fixed robust control during identification 

phases is closer to the simultaneous identification and control 

objective (and may avoid the stall mechanism noted in the previous 

section) but does spread the W_. spectra again severely limiting 

the benefits of the LPF. Another seemingly applicable concept is 

28 

that of adaptive orthogonal filtering . This idea is incorporated 
by appending to (8) and (13) additional uncontrolled "modes" 
intended to absorb the spillover effects* These modes would 
require time-varying dynamic descriptions meaning that the p and 
p used in estimating their difference equation parameters would need 
to be significantly larger than those used in estimating the time- 
invariant modal difference equation parameters* The compensatory 
ability of these additional model modes seems limited due to the . 
assumptions necessary in adaptive orthogonal filter development. 
Therefore none of these suggestions appears wholly satisfactory. 

The conclusion is that currently the observation spillover problem 
remains unsolved yet requires resolution for broad applicability 
of the proposed adaptive modal control scheme for BPS. 

Conclusion 

This paper begins with revision of a previously originated 

4 

strategy for adaptive modal control of DPS and concludes with the 
confrontation of the spillover problem, which is extremely severe. 

In support of the simulation evidence provided it has been proven 

5 

elsewhere that in the absence of observation spillover and with the 
use of the eigenvector 'expansion, adaptive* modal control of DPS 
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is as viable as lumped parameter system indirect adaptive control- 

g 

It can also be shown why a stable simultaneous identification 
and control scheme similar to that imbedded in the AMCD example 
fails in the presence of observation spillover or nonorthogonal 
expansion. With the necessity of a reduced order model (N > M) 
the goal of globally stable adaptive DPS controller convergence 
appears too stringent. Work is in progress to relate observation 
spillover bounds to parameter identification bounds. Such efforts 
are directed toward delineating the detrimental influence of observation 
spillover and the possibility of allowable behavior despite its 
presence rather than toward its removal. 

29 

Suggestions have also been forwarded to remedy the difficulty 

of selection. For general spacecraft, the fundamental assumption 

of eigenvector availability for (1) is overly optimistic. Slight 

inhomogeneities in the AMCD can lead to significant coupling of 

30 

the Fourier expansion "modes" . Such coupling, if incorporated 
in the system model, disallows the parallel computation objective 
for the identifier and "modal" controller. The specter of the 
necessity to recursively estimate both the basis functions and their 
associated dynamics raises questions far beyond the current scope 
of these efforts. However, these are issues that ultimately must 
be addressed. 

The development of an adaptive controller applicable to DPS 
requires examination of both indirect and direct adaptive control 
concepts in a necessarily reduced-order model format. Both approaches 
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are susceptible to spillover degradation. A judicious mixture of 
robust control, gain-scheduling, on-line versus off-line identification, 
specific "optimal ,, objectives versus simpler damping requirements, 
and local versus global convergence behavior will be required 
in subsequent efforts. Further research, as in any emerging field, 
will better identify the weaknesses and strengths' of proposed 
approaches to adaptive control of DPS and uncover additional concerns 
requiring further original developments. 
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On-line adaptive control is essential for Large Space Structures 
(LSS) where the modal parameters are poorly known, due to modeling 
error, or changing, due to variable configurations. It is especially 
important that such adaptive controllers produce stabilizing controls 
during adaptation due to the small damping present in LSS. However, 
any such controller must be based on a reduced-order model of the LSS. 
The spillover from the unmodelled residuals, as well as the modeling . 
error, can- deteriorate the performance of the adaptive controller and, 
if uncompensated, this spillover can defeat the whole purpose of the 
adaptive control. 

This paper investigates adaptive control for LSS using direct and 
indirect schemes and points out the mechanisms whereby observation 
spillover can upset the stability of the controller. The framework for 
nonadaptive control of LSS is reviewed and many of the generic problems 
of adaptive LSS control are pointed out within this framework. These 
generic problems must be overcome for successful operation of adaptive 
LSS control . 


1.0 INTRODUCTION 


This paper deals with the basic problems inherent in adaptive con- 
trol of large space structures (LSS), such as satellites and spacecraft, 
where the structural parameters are poorly known or slowly time-varying. 

With the advent of the Space Shuttle Transportation .System, it has 
become possible to conceive very large spacecraft and satellites which 
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could be carried i-nto space and deployed, assembled, or manufactured 
there. Such LSS would serve a variety of civilian and military needs 
[ 13 , [2], including electrical energy generation from the solar power 
satellite - a structure nearly the size of Manhattan Island - to be 
constructed in space and operated in earth geosynchronous orbit [3], The 
control technology- needs for such LSS have been discussed in a variety 
of articles, e.g., [4], [5], and the developing LSS control theory and 
technology has been surveyed, for example, in [6]-[10]. 

The size of these structures, their low rigidity, and the small 
damping available in lightweight construction materials combine to make 
LSS extremely mechanically flexible. In theory, LSS are distributed 
parameter systems whose dimension is infinite; however, in practice, 
their dynamics are usually modeled by Targe scale systems based on 
approximate elastic mode data. Active control schemes for LSS are often 
required to meet stringent requirements for- their shape, orientation, 
alignment, and pointing accuracy. Such active control is limited by the 
capacity of the, on-board control computer, the modeling inaccuracy in cur- 
rent finite element computer codes for analyzing structural dynamics, and 
available control devices (actuators 'and sensors); therefore, the. con- 
troller must be based on some reduced-order model (ROM) of the LSS. 
Fundamental problems of LSS control include; 

(1) selection of appropriate modes to control for des.ired system- per- 
formance; 

(2) development of ROM for analysis and controller design; 

(3) computation of system model and control parameters; 

(4) controller design with multiple distributed actuators and sensors; 

(5) the number and location of sensors and ‘actuators for efficient . 
control-; 

(6) the- effect of, and compensation for, residual ’ (unmodel led) modes 
and modeling error on the- closed-loop system .performance; 

(7) adaptive, and self-tuning controllers for LSS with poorly known or 
changing parameters and configurations. • 

Item (7) is the basic topic of this paper but it must be considered in 
the context of the other items with which it is completely intertwined. 

The need for adaptive control in LSS arises because of ignorance of ' 
the system and changing control regimes . The former occurs as 
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(a) ignorance of the system structure and order, and (b) ignorance of the 
■system parameters; the latter occurs because of changing configuration 
of the- LSS. Changes in configuration may be due to construction in-space, 
thermal distortion, or reori entation of subsystems, e.g., rotating solar 
panels or sunshields; these changes usually produce slowly time-varying 
parameters. Ignorance of the LSS system structure and order -is due to the 
fundamental problem of modeling a distributed parameter system,, e.g., 
faulty physics, reduced-order models, and ignored nonlinearities; this 
means that the order of the ROM is lower than that of the actual LSS. 
Ignorance of the system parameters , while directly related to the system 
structure, is due to the inherent modeling error present in even the best 
structural analysis computer codes and to the limitation of testing huge, 
lightweight LSS on earth; this produces constant but poorly known system 
parameters. There is a very clear need for an adaptive LSS control metho- 
dology that' can begin with the best available computed parameters and 
self-tune its way toward the correct parameters while stablely controlling 
the LSS and, .possibly, adapting to variable configurations. 

A great variety of adaptive control schemes exists for lumped param- 
eter, small scale systems [11]; in particular, model reference adaptive- 
methods have achieved a great amount of success in producing stable, 
convergent adaptive controllers, and adaptive observers for systems whose 
structure is known and whose parameters are constant but poorly known 
or slowly time-varying, e.g., [12]-[25]. Adaptive schemes may be direct , 
i.e., the available control parameters are directly adjusted (adapted)- to 
improve the overall system performance, e.g., [25]- [26] , or indirect , i.e., 
the system parameters are identified (based on the assumed system struc- 
ture) and the control commands are generated from these parameter esti- 
mates as though they were the actual values, e.g., [20], [24], [27]. 

The abundance of adaptive control methods is overwhelming and an 
understanding of the interrelationships and structural commonality of these 
methods is desperately needed; see, e.g., [28], [29], for some beginnings 
in this direction. Furthermore, • the use of such methods on distributed 
parameter or large scale systems, like LSS, is greatly limited by the ROM 
problem - the. adaptive scheme must be based on a ROM of the actual, system 
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and, hence, the order of the model is, and must remain, substantially 
lower than the controlled system. In addition, it inessential that the 
LSS adaptive controller provide a stabilizing control in such highly 
oscillatory systems as LSS . 

This paper develops a framework for LSS adaptive control problems and 
points out generic problems in the use of the most natural direct and 
indirect adaptive approaches. In other forms, these problems will haunt 
every use of adaptive control on LSS and must be solved before the valu- 
‘ able benefits of adaptive control can meet the needs of this new applica- 
tion area . A few preliminary attempts at adaptive control for specific 
distributed parameter systems or LSS have been made in [30]~[36]; also, 
for the corresponding parameter identification problem for distributed 
parameter systems, see [37]. 

2.0 NONADAPTIVE LSS CONTROL 

Following [6j, the LSS may be described by the partial differential 
equation: 

m(x) u tt (x,t) -s- D q u t (x,t) + A q u{x,t) = F(x,t) (2.1) 

where u(x,t) represents (possibly, a vector of generalized) displacements 
of the structure ft off its equilibrium position due to transient dis- 
turbances and the applied force distribution F(x,t). The mass distribu- 
tion m(x) is positive and bounded on ft. The internal restoring forces of 
the structure are represented by Aq u where Aq is an appropriate differen- 
tial operator with domain D(A Q ) defined in a Hilbert space Hq with 
inner product (.,.)q. In most pro (foo t cases, Aq has discrete spectrum, 
i . e . , 

A 0 * k ■ <4 \ (2-2) 

where ou^ are the mode frequencies of vibration and (^(x) are the mode 
shapes . The damping term Dq u t is generated by an appropriate A Q -b6unded 
differential operator and may represent gyroscopic damping as well as the 
very small (^a critical) natural damping expected in the LSS. 

The applied force distribution is given by 

F(x,t) = F c (x,t) + F 0 (x,t) 


(2.3) 
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where Fg represents external disturbances and Fg represents the control 
forces due to M actuators: 


• M 

F c Cx,t) = B 0 f = X 
1-1 


(2.4) 


where are the actuator influence functions (usually point devices) and 
f.- are the control commands. Observations are produced by P sensors: 

y = C 0 u + C 0 u t (2 ‘ 5} 

I 

where y.(t) = (c.,u) n + (c.,u.) n for 1 < j < P with c. being the position 
sensor influence functions and c( the velocity sensor ones (usually point 

V * 

devices). 

The state variable form of (2.1) and (2.3)-(2.5) is obtained by 
taking 

T 

v(x,t) = [u(x.t-), u t (x,t)] 
in H = 0(Ag) x Hg with energy norm: 

|| v || 2 « (mu t , u t ) + (Ag u, Ag u) (2.6) 

This produces 


v. = A v + Bf; v(0l * w g 


(2.7) 


y = Cv 

where we consider the disturbance-free case (Fg = 0) and define B = [0 Bg]" 1 ", 

~0 I~ 

his distributed parameter system is 


C - [C 0 Cg] and A = 


L‘ a o 


-D, 


very oscillatory in the sense that the semigroup U(t) generated by A has 
very 1 ittle damping: 


U(t) II 1 M n e " £t for t > 0 


( 2 . 8 ) 


where e >_ 0 and small and Mg >_ 1 . 

The desired performance of the actively controlled LSS greatly effects 
the design of the controller. Many desirable properties of the active 
structure can be obtained with constant feedback gains applied to the sys- 
tem state v(x,t); such solutions arise for regulator problems and sta- 
bilization (pole placement) problems for LSS. However, the full (infinite 
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dimensional) state v is never available from a distributed parameter sys- 
tem; only the P sensor outputs y are available. 

Impl ementable controllers for LSS (and most distributed parameter 
systems) must be based on finite dimensional on-board control computers 
which process the sensor outputs y and produce control commands f; thus, 
a reduced-order model (ROM) of the system (2.7) must be used for the 
controller design. A ROM can be obtained by projecting the system (2.7) 
in H onto an appropriate finite dimensional subspace the projections 
P (onto H^) and Q (onto the residual subspace ) are usually, but not 


always, orthogonal. Let v^ = Pv and v 


“ A. 
N i 


V N + A NR V R + B N f 


V R = A RN V N + A R V R + E R f 


V N + C R V R 


Qv and, from (2.7), we obtain: 

(2.9) 

(2.10) 

( 2 . 11 ) 


y = C 

where = PAP, A 

called control and observation spillover ; the terms A 
model error . The ROM for this system is given by (2.9) and (2.11) with 
A NR = 0 and C R = 0: 


^NR = BA< ^’ = PBs e ^ C * ”^ e terms ^ and ^R V R are 

'NR V R ar8 ca ^ ed 


V N “ A N V N + S N f 

y = C N V = % 


( 2 . 12 ) 


The ROM state v^ and the residual state v^ form the true system state v 
with total energy |jv|| 2 given by: 

V = V N + V R 

N V H 2 = l|v N ll 2 + 1! v R |j 2 (if projection is orthogonal) 

All implementabl e controller designs based on any ROM must be evaluated 
in closed-loop with the actual LSS (2.7),and it is in this evaluation that 
the effects of model error and spillover due to the residuals become 
apparent. 

If the actual mode shapes are known, the modal ROM is a sensible 
choice : 


~{sp <j>j| ... tp^} 

and the model error terms A^ and A^ become zero. Of course, any collec- 
tion of modes could be used; usually, the most easily excited or critical 
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ones will be chosen. However, in many cases, the partial differential 
operator A is too complex to provide closed-form mode shapes. Instead 
finite element approximations of the mode shapes are computed (e.g., 
via NASTRAN) / and these approximate mode shapes can be used to form the 
ROM for controller design; note that some model error is present when 
these approximations are used. Henceforth, we will a ssume the actual 
mode shapes are available to simplify the discussion but much of our 
analysis remains valid for approximate mode shapes and other types of 
ROM. 


Modern modal control (MMC) for LSS, as developed in [38], uses the 
modaT (or approximate modal) ROM and develops a controller consisting 
of a state estimator based on the ROM and a constant gain control law: 


? N = A N ? H * 8 N f + 
y ~ V N’ V N^ = ® 


(2.13) 


and 



(2.14) 


This controller design requires the ROM (A^, B^, CjJ to be controllable 
and observable for the calculation of control and estimator gains G^, 
K^ t . These conditions, in modal terms, provide insight into the number 
and location ^f actuators and sensors. From [38], (A^, B^, C^) is con- 
trollable and observable, when position sensors are used (Cg = 0), if 
and only if 


(1) min(P,M) >_ max .mode frequency multiplicity in the ROM 

(2) each sub-block of 8^ and associated with a mode frequency of 
multiplicity must have rank at least equal to cs^,. 

Similar results hold for other types of sensors, e.g., velocity, acceler- 
ation, or mixtures of types [38], [39]. These results are easy to in- 
terpret in terms of the mode shapes, e.g., if no repeated frequencies 
exist, then the above result says that a single actuator and sensor, 
not necessarily collocated, will do the job as long as neither is located 
at any of the ROM mode shape zeros. Since LSS have many symmetries and 
rigid body modes, it is not often that a LSS control problem will have a 
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controllable observable ROM with only one pair of devices; this has con- 
sequences for the adaptive control problem to be discussed later. 

Let the estimator error -e^ = v^ - v^ be defined and, from eqs. (2.9) 
-(2.11) and (2.13)-(2 .14) , obtain 


*N " ^ A N + B N V V N + B N G N e N 
= (A N ' K N V e N + K N C R V R 


(2.15) 

(2.16) 
(2.17) 


This shows the effect of spillover on the closed-loop system: even 

though G^, A^ - C^, and A^ are stable, the closed-loop system 

need not be stable. When either control or observation spillover is 
absent (B^ = 0 or C R = 0), then stability is assured; otherwise, spill- 
over causes pole-shifting and can induce instabilities [38], [40]. 

Bounds on the destabilizing effect of observation and control spillover 
were produced in [38] 'and can be extended to the case where some model 
error and small nonlinearities are present [41]. Such bounds give an 
indication of how much spillover the closed-loop system can tolerate. 

A variety of methods have been suggested to reduce spillover [6]. 

One obvious way would be to prefilter the sensor outputs, with a bandpass 
filter, to substantially reduce observation spillover. This alleviates 
the worst pole-shifting problem; bounds on system performance with' con- 
trol spillover alone can be found in [42]. Note that the post-filter 
of the controller outputs could do the same job by reducing control 
spillover; this interchange of filter and controller is possible due to 
linearity and time-invariance. The trade-off with this means of reducing 
spillover is that the prefilter introduces phase distortion which can 
have a destabilizing effect of its own. Therefore', a very high order 
filter may be required to keep the phase distortion acceptable; phase- 
locked-loop quadrature filters may be another solution [43]. E ven in 
nonadaptive LSS control, the spillover and model error problem is a funda- 
mental one . 

Finally, we should note in this section that digitally implemented 
controls would be based on discrete-time versions of the distributed 
parameter system (2.7). One such version is obtained by using a uniform 
time step At: 
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v(k + 1 ) = $v(k) + E 
y( k) = Cv(k) 


(2.18) 


At 

where $ = U(At) and E = E Q U(t) dx B and the control command is 


a constant fj, over the interval (k - 1 ) At <_ t < k At. Other versions of 
this could be obtained with nonuniform time steps. When the ROM procedure 
of projecting onto the subspace is used, we obtain the discrete-time 
versions of (2.9)-(2.11):‘ 


V k ' 

+ D 

= * N v N (k) + 

^NR + E N T k 

(2.19) 

v R (k 

+ 1) 

= $ RN V k * 


( v R (k > + E N f k 

(2.20) 


y(k) 

= C N v N< k > + 

C R 

v R (k) 

(2.21) 

H n 1 = 

; the 

modal subspace. 

the above become: 


v N (k 

+ 1) 

= v N eo + 

E o 

6 N f k 

(2.22) 

v R (k 

+ 1) 

= $ R v R (k) + 

E 0 

B R f k 

(2.23) 


y(k) 

= C, v N (k) + 

C R 

v R (k) 

(2.24) 


where (2.22) is the same as that obtained by directly discretizing the 
ROM in (2.12); if the exact mode shapes are not available, these two 
discretizations may yield different results. In addition, the sampling 
process can alias residual modes and increase observation spillover and 
the .zero-order hold process can spread-out the control command signal 
spectrum and, hence, increase control spillover by increasing the energy 
in the residual mode spectrum; this has been observed and investigated 
in [44], Therefore, the time discretization is a very important factor 
in the design of implementable LSS controllers . , 


3.0 TOWARD ADAPTIVE CONTROL OF LSS 


In order to design MKC, or other controllers, for LSS as proposed in 
the previous section, we must have knowledge of the ROM parameters (A^, 
B^, C^). These parameters are obtained from modal data; they are the 
mode frequencies for A^ and the mode shapes at actuator and sensor loca- 
tions for and C^,, respectively. This data is required for three 


reasons : 



TO 


(1) to determine controllability and observability of the ROM and, 
hence, to help locate control devices effectively; 

(2) to design control and estimator gains; 

(3) to use in the state estimator's internal model. 

However, we have noted in Sec'. 1.0 the sources of error for this data; 
consequently, a need arises for an adaptive version of the MMC of Sec. 

2 . 0 . 

The most logical and reasonable procedure to obtain adaptive con- 
trollers for a LSS seems to be the following: 

Procedure for Adaptive LSS Control 

(a) choose a "nice" reduced-order model (ROM); a modal ROM 
would be the obvious choice; 

. (b) use your "favorite" lumped parameter adaptive control 
scheme; 

(c) design the adaptive controller as though the ROM were 
the actual LSS to be controlled, i.e., ignore the un- 
modeled residual part of the structure; 

(d) use this adaptive controller in closed-loop with the 
actual LSS and hope for the best. 

There is nothing wrong with Tollowing this procedure as a best first guess 
- in a way, there is little else that one can do to produce an implement- 
able adaptive LSS controller. 

In some cases, spillover is sufficiently small or enough other mathe- 
matical structure is present in the system, e.g., a high level of damping 
in the distributed parameter system [31] or low level of performance re- 
quired from the controller (increased damping via direct velocity feedback) 
[32], to allow the adaptive controller to operate. However, these situa- 
tions are rare with LSS and one would not like to count on the "generosity 
of nature" or the temporary suspension of Murphy's Law as part of the 
above design procedure. Consequently, we would add the following items to 
that procedure: 

(e) analyze computer simulations of higher-order models of the 
LSS in closed-loop with the adaptive controller based on 
the lower-order ROM (e.g., simulate more modes than you 
plan to control ) ; 
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(f) investigate the specific mathematical mechanisms whereby 
the residual (unmodeled) part of the LSS couples into a 
given adaptive control scheme (e.g., find out where and how 
spillover affects the adaptive controller)'; 

(g) obtain mathematical results on the amount of spillover 
and/or model error that can be tolerated in the closed- 

loop system and still provide adequate adaptation and control; 

(h) develop spillover and model error compensation schemes to 
augment the adaptive controller when the residuals cannot 
be tolerated (e.g., when the conditions of (g) are not sat- 
isfied); 

(i) recheck (g) with this compensation in the closed-loop system. 

We believe that, within the basic framework of LSS control as de- 
scribed in Sec. 2.0, this Augmented Procedure (a)-(i) will go a long way 
toward revealing the problems of adaptive LSS control (and indeed, most 
other adaptive control situations where control must be based on a ROM) 
and will help to- focus needed attention on these crucial issues. For 
example, although Step (e) would be done most ljj££l^-ajL some point in 
the system development phase of any project as^tne construction and 
operation of a LSS, often it is done much too late and the design is 
"set in concrete (or in this case, graphite-epoxy)"; the other steps 
(f)-(i) may not be done at all. Yet, ignoring the effects of the resid- 
ual unmodeled LSS can produce some disasterous behavior in the adaptive 
controller; this was pointed out quite clearly in the LSS example in [36], 

When an adaptive scheme is applied to such a lightly damped, oscil- 
latory system as a LSS, the stability of the closed-loop system during 
adaptation is a necessity ; therefore, we do not view convergence. and 
(global) stability results as luxuries for adaptive LSS control and shall 
only consider appropriate those lumped parameter adaptive schemes for 
which such results, are available. However, even a globally stable adaptive 
scheme may prove to be unstable when it is used in closed-loop with the 
actual LSS instead of the ROM on which it was based . This is not a failure 
of the, adaptive scheme; it is a failure to satisfy the mathematical hypoth- 
esis of the stability result associated with the scheme. 
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In the rest of this section, we shall study the use of two very 
general adaptive schemes which seem to illustrate the problems and poten- 
tial of adaptive LSS control: the indirect schemes of [23]-[24], which 

use an ad aptive observer and operate in continuous time, and the discrete- 
time, dir ect or indirec t , schemes which are based on an Autoregressive 
Moving Average ( ARMA) Model of the controlled plant, e.g ., [25J-[27]. 

These approaches represent a good cross-section of available lumped 
parameter adaptive control schemes which have been shown to possess the 
desired stability properties. We emphasize that the. point of this secr 
tion is not to criticize or slander these schemes; rather, we mean to 
point out where the hypotheses of their stability results are violated, 
and must be modified, when we attempt to use them on LSS. We feel that 
consideration of these approaches within the context of the Augmented 
Procedure for Adaptive LSS Control (a)-(i) will illustrate the generic 
difficulties in the application of existing, well-behaved, lumped param- 
eter adaptive control schemes to LSS. 

3.1 Multivariable Systems Converted to Scalar Systems 

The results of many stable adaptive schemes, e.g., C'13]-[14], [23]-[24], 
are limited to a single actuator and/or a single sensor; yet, we have seen 
in Sec. 2.0 that most LSS co:.tr.ol problems will involve multiple actuators 
and sensors. One way to deal with this (although, admittedly it has its 
drawbacks) is to convert the controllable observable LSS problem via out- 
put feedback int.o one that is controllable and observable from a single 
actuator and/or single sensor; this can be done with almost any output 
feedback gains [45]-[47]. These gains would have to be based on the best 
available calculated ROM data and the designer must hope that they wi'll 
continue to do their job during adaptation. 

The output feedback modifies the original system (2.9)-(2.11) to be- 
come: 

v N = (A n + B, H n C n ) v n + (A nr + B n H n C r ) v R + b N f 

*R = ( A RN + B R - H N V V N + ( A R * B R H N C R' V R + b R f 
y = C N V N t C R V R 

where is the output feedback gain matrix, b^, b R , c^, c R are vectors. 


(3.1) 

(3.2) 

(3.3) 
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and f,y have been renamed. Let A^ + be A^, etc. and we have 
that^new ROM (A^, b N> c^) is a controllable, observable single input, 
single output system and (3.1) and (3.2) become: 


*N = A U v N + A NR v R + b « f (3 ' 4) 

V R ” fl RN V N + A R V R + b R f * 3-5 ' 


3.2 Indirect Adaptive Controller Design 


We apply the design of the adaptive controller in [24] directly to 
the ROM consisting of (.3.4) and (3.3) with the assumption, for now, that 

A^ - 0 and c^ - 0 . 

The control law is given by 

f(t) = gj v N (t) + f.(t) (3.6) 

where g^ is a constant gain vector, f is a "sufficiently rich" external 
signal (more about this later), and v^ is derived from the following • 
adaptive observer (or state estimator): 


A 

v. 


= F v„ 


gy 


hf; 


V N (°) • v 0 


(3.7) 


where F is an arbitrary, stable matrix and g, h are unknown parameter 
vectors. The appropriate m atching cond i tio ns are: 


F + G*cJ = A n 


h* « b 


v 0 = V N (0) 


$m) 

(3.8) 


where g*, h* are constant. 

Let p* = [g*^ h* T v n (0) T j T and note that 

p* = F p*; p*(0) = p Q (3.9) 

where F = diag [0 0 F] and, when g = g*, h = h*, -we have from (3.8) and 


v N (t) = [M(t) I N ] p*(t) (3.10) 

where M(t) e [I N y(x) 1^ f(x)] dx and we have used the fact 

that h, f are scalars . 
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Mow (3.7) can be rewritten as 


v N (t) = CM(t) I N ] p(t) 

M(t) = F M(t) + [I N y(t) I N f(t)] 
M(0) = 0 


(3.11) 

* 

(3.12) 

(3.13) 


where M(t) is as defined in (S.lO^and we have yet to produce the adaptive 
law to generate p(t). This adaptive law is given by the following: 


p(t) = F p(t) - a(t) <"R(t) p(t) + r(t)] 
where p(0) = p Q is arbitrary and 


(3.14) 


'R(t) = -q R(t) - F T R(t) - R(t) F + [M(t) IJ T C N cj [H(t) I w ] 


and 


R( 0) = 0 


r(t) = -q r(t) - F T r(t) - [M(t) I N ] T c M y(t) 


N- 
(3.15) 


(3.16) 

r(0) = 0 

where a(t) is the adaptive gain and the constant q is chosen to exceed 
twice the absolute value of the real parts of the eigenvalues of F. The 
adaptive gain is chosen so that 

a(t) ■ Y * llitt)! I (3.17) 


where 

p(t) = -X y(t) + (N/2) J ‘ (iy{t)| + |f(t)|) (3.18) 

with X positive and F + < 2X 1^. 


3.3 Convergence Results: What Goes Wrong? 

All of the above is exactly as stated in [24] where it is also shown 
in Appendix I and II that 

y(t) > | | M(t) | j (3.19) 

and, with f(t) sufficiently rich in frequencies, there is a t^ such that 
R(t) >_ p > 0 for all t >_ tj (3.20) 
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In addition, it is shown in Appendix III that 

R(t) p*(t) + r(t) = 0 (3.21) 


This is very crucial to the stability results of [24] and it is here that 
observation spillover (i.e., the fact that we are using a ROM) appears - 
(3.21) is not valid when c R f 0; however, 

R(t) p*(t) + r(t) - A R (t) (3.22) 

t ~T 

where A R (t) i J e”^ q + F ) - T ) [m( t ) i n ] T c n c R v r (t) dr 

Let e N (t) = v N (t) - v N (t) and Ap(t) = p(t) - p*(t); then 

e N (t) = [M(t) I M ] Ap(t) (3.23) 

and 


Ap(t) = [F - a(t) R(t)] Ap(t) - a(t) A R (t) (3.24) 

Consider V(t) = Ap(t)"*" Ap(t) and we obtain: 

V(t) < -y P V ( t ) - 2y A R (t) Ap(t) (3.25) 

where V(0) = ApJ Ap Q and Ap Q = p Q - p*. 

This follows [24] except for the additional term in (3.25); also 
note that 


e N (t) ! | 2 < V(t) [1 + u(x)] 


(3.26) 


Let Av N (t> = v N (t) - v*(t) and Av R (t) = v R (t) - v*(t) where v^ and v* 
represent the idea*l states of (3.4) — (3.5) when the parameters are 
exactly known: 


V H “ ^ A N + b N 9 N } V N + A NR V R + b N f c 

V R = <*RN + b M * V N + A R V R + b R f c 
with v*(0) = v N (0) and v*(0) = v R (0) . 


(3.27) 

(3.28) 


When the implementable adaptive control law (3.6) is used, we obtain: 

(3.29) 


AV N ' ( A N + b N g N^ Av N + A NR AV R + b N 9 I e N 




AV R = ( A RN + b N V AV R * A R AV R + b R 9 N e 


N 


(3.30) 


we have the following result: 
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THEOREM 3.1 : Assume 

(1) c R = 0 

(2) A N + b N gj stable 

(3) A R stable 

(4) A nr =0 

(5) f (t) sufficiently rich (i.e., it has at least 3N/2 
distinct frequencies) 

O 

Then there is a S > 0 such that for all j |aPq| | ± 6: 

(a) e^(t) ls bounded and (eventually) vanishes with an 
arbitrary exponential rate 

(b) 1 im Av N (t) = 0 

(c) lim Av D (t) = 0. 
t-*» K 

Therefore, even though the closed-loop system with the adaptive 

controller is highly nonlinear, it is stable while the adaptation 

is taking place. In particular, (1), (3), and (4) are satisfied 

if there is no observation spillover (C R = 0) and no model error 
(A^ R = 0) and some damping in -the residuals (i.e., A R is stable). 

Proof: This result follows from the results of [24] because A n (t) = 0 in 

K 

(3.22) when c R = 0. The stability of (3.29)-(3.30) is determined by 
A n + b N gj and A R when A NR = 0. Also, if C R = 0, then A NR = A^ R , A R = A R , 

and c R = 0 since it is a row vector of C R . # 

Note that the stability of the system with adaptive control is 

determined by that of (3.27)-(3.28) - the ideal case where the parameters 

are known and the external sianal f is applied. This is natural, since 
— — c 

adaptation cannot take place without f present; however, after adapta- 

’ V 

tion, we would most likely want to turn off f c . In addition, we could 

choose the f signal so as not to excite the residual frequencies whenever 

sufficient spectral separation is present. Still, we would need to turn 

on f„ now and then, in order to "tune-up" the controller, 
c 

3.4 Spillover Compensation for the Indirect Adaptive Controller 

The above result is merely a slight extension of the results of [24] 
to a special case of the adaptive controller based on an ROM instead of 
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the full system. However, i-t does suggest that some form- of compensation 
should be used to eliminate the observation spillover. Such compensation 
must be essentially independent of the parameters of the ROM; most 
methods of spillover reduction require- knowl edge of the ROM (and some 
residual) parameters. 

One approach to spillover compensation already suggested in Sec. 2.0 
is prefiltering the sensor outputs to remove or greatly reduce the term 
v^(t). Such prefiltering can be achieved with low-pass or band-pass 
filters when the ROM frequencies are known and separated from the residual 
frequencies. However, the modal frequency data is part of the poorly 
known parameter data. 

In an attempt to resolve this predicament, we could try using phase- 
locked loop (PLL) based filters with the center frequency of each loop 
tuned to the best approximation available of the corresponding ROM The 
PLL will adapt itself until' it tracks the actual mode frequency and, 
after "lock-on," it will behave as a narrow band-pass^l inear filter which 
tunes out the observation spillover from the other frequencies [43]. Of 
course, sufficient spectral separation must be present, the calculated 
values of the ROM mocal frequencies must be sufficiently good, the dis- 
tortion introduced by the filter must be sufficiently small, and the 
adaptive controller rust not shift the poles around too much. Thus, the 
PLL prefilter is’not a panacea/ But, it might work to reduce spillover 
and, if it does, it has the added advantage that its output could also 
reveal better estimates of the modal frequencies; this would take some 
of the load off the adaptive observer. If the adaptation mechanism causes 
too much pole shifting, the ROM -frequencies could be excited via f and 
identified in' open-lcop by the PLL filters before the adaptive controller 
is turned on. Note that modal frequency data is usually better known, 
via computer approximation, than modal shape data; hence, this approach 
might not be unreasonable. Another possible, but untried, approach to 
spillover compensation might be an adaptive version of the orthogonal 
filter, in [48]. 

Note that some prefiltering (and - postfiltering) always takes place 
due to the bandwidth limitations of the sensors (and actuators).' Whether 
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A problem that arises with the use of the results in [24] and their 

T 

modification to LSS is that the constant feedback gains must be cal- 
culated, in advance, to stabilize + b^ gjjj. It would be better if 
these gains were adapted along with the parameters in the observer. Of 
course, after adaptation has taken place, they could be recalculated 
from the "tuned-up" parameters, but, in some cases, the adaptation phase 
is never over, e.g., slowly-varying parameters. Other approaches could 
be used for adaptive pole-placement, e.g., [20], [21], but these also 
have their limitations. 


3.5 ARMA-Gettin 1 

Many discrete-time adaptive control schemes depend on an Auto- 
Regressive Moving Average (ARMA) representation of the plant in discrete- 
time, e.g., [25]-[27]: 

N N 

y(k + N) “ 2? a r y( k + r - 1) + 2. S r f ( k + r ' (3.31) 

r-1 r=l 

for some N and appropriate matrices a f , 6 r . What the ARMA says is that, 
after N time steps, the present output is only related to the past N out- 
puts and inputs. Existence of an ARMA is directly related to the finite 
dimensionality of the plant (N is usually that dimension) and is obtained 
using the Cayley-Hamilton theorem for matrices. For LSS, only a " quasi- 
ARMA" can exist; these were considered in detail in [49]. From the 
Appendix of [49], we obtain the quasi-ARMA for the LSS (2.18) or (2.19)- 
( 2 . 21 ): 

N N 

y(k + H) * 2 <* r y(k + r - 1) + £r r E^ f(k + r - 1) + R(k) (3.32) 
r-1 - r=l 


N 

R(k) = C R v R (k + N) + 2! v R ( k + r - 1) 

r-1 


A = 
r 


^NR 


“r C R 


where 
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and T is defined in the Appendix of [49]. Since R(k) = 0 when C D = 0 
and A r = 0, we have the following result: 

THEOREM 3.2: When the observation spillover (C^) and the 

model error term are both zero, the quasi-ARMA (3.32) 

.is a true ARMA for the LSS (2.18) / and any stable adaptive 
scheme based on this ARMA will be globally stable when used 
in closed-loop with the actual LSS (2.18). 

When the rather stringent hypothesis of Theo. 3.2 is not satisfied 
(as it may not be in practice), any adaptive LSS control scheme based on 
the quasi-ARMA (3.31) must ignore R(k) in order to be imp! ementabl e . How- 
ever, R(k) is the term where the residual effects - spillover and model 
error - enter the scheme and can cause instability. Again, as in Sec. 

3.4, prefiltering or other compensation might be tried in an attempt to 
reduce or eliminate this term. 

4.0 CONCLUSIONS 

In an attempt to point out the crucial issues and generic problems 
associated with adaptive control of large aerospace structures (LSS), we 
have reviewed the framework for nonadaptive LSS control (Sec. 2.0) and, 
within this framework, have proposed a general procedure, based on 
reduced-order models (ROM) of the LSS, for obtaining and assessing the 
problems of adaptive LSS controllers (3.0). In addition, we have taken 
a closer look at the use of certain well-known, lumped parameter, stable 
adaptive- control schemes in this procedure. Taking these schemes as 
representative of the basic ideas present in all lumped parameter adap- 
tive control approaches, we have obtained corresponding LSS adaptive 
controllers and found the following generic problems associated with 
adaptive LSS control : 

(1) LSS are distributed parameter or large scale systems; 
therefore, the plant dimension is always larger than 
the dimension of the adaptive controller, which must be 
based on a ROM; 
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(2) LSS control must often be done with more than one actua- 
tor and sensor; conversion of multivariable to scalar 
systems via output feedback introduces problems of its 
own (e.g., stability of the residuals); 

(3) LSS control problems are often- non-minimum phase due to 
nonco.llocated actuators and sensors; 

(4) Interaction of the residuals with the adaptive control- 
ler may negate the stabilizing properties of the con- 
troller due to observation spillover; this interaction 
is much worse due to the nonlinear nature of adaptive 
control ; 

(5) Methods of spillover compensation for LSS often require 
knowledge of the ROM (and some residual) parameters - 
t.he very data that are poorly known; 

(6-) The adaptation mechanism may shift the closed-loop fre- 
quencies around; this counteracts the benefits of any 
prefilt-ering unless sufficient spectral separation is 
maintained; 

(7) Indirect adaptive controllers need sufficient excitation 
from an external signal f ; however, this signal may 
substantially excite the residuals. 

(8) Discrete-time adaptive controllers can only be based 
on quasi -ARMA rather than strict ARMA representations 
of the LSS; this may negate the stability properties of 
such a controller. 

Stable adaptation is essential for such highly oscillatory systems 
as LSS, yet our preliminary stability results, Theos . 3.1 and 3.2, both 
require that observation spillover be somehow completely eliminated 
before it reaches the adaptive control logic; certainly, this is not 
an easy thing to do in general! Perhaps, global stability is t,oo much 
to ask for LSS adaptive control because observation spillover will 
always be present to'some degree in LSS control. However, it seems rea- 
sonable to hope for the development of spill over bounds to give some 
idea of regions of stability for the successful operation of LSS adaptive 
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control. Some comparison should be made between stable adaptive con- 
trollers based on ROM and stable, robust control schemes, e.g., [50],, 
[-51]. . _ ' 

In closing, we would like to say that it is not our - intent to present. 

. i ■ i r^i->r-»rTrTi -■ Mri — — • •• *-' r 

a- gloomy picture for the application of adaptive control to LSS. In 

— . . . . : — ' — TT i nr . ir .itjiey 

fact, the need for adaptive control in LSS is already becoming quite 
clear, and recognition of this need comes, for a change, at an appropri- 
ate time - before any LSS have been built and put into space. However, 
the development of. adaptive control for LSS will not take place overnight 
and will not be done by one or two people. Consequently, what we have 
tried to stress here for the interested researcher are some of the funda- 
mental problems that arise and the basic steps which need to be taken 
toward the goal of successful adaptive control of LSS. In the long run, 
we have high hopes for the success of this endeavor and -we expect' that 
adaptive control theory will profit by its association with large aero- 
space-structures, as well.. 
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